PROCEEDINGS 

OF SCIENCE 



Compatibility of global NPDF analyses of neutrino 
DIS and charged-lepton DIS data 



o 

(N 

> 
O 



Or 
43 



> 
in 



X 



K. Kovarik 

Karlsruhe, Institute of Technology, D-76128, Germany 



E-mail: kovarik@particle . uni-karlsruhe . de 



I. Schienbein 

LPSC, Universite Joseph Fourier Grenoble 1, CNRS/IN2P3, INPG, Grenoble, F-38026, France 
E-mail: 



schien@lpsc.in2p3.fr 



F.I. Olness 

Southern Methodist University, Dallas, TX75275, USA 



E-mail: olness@physics . smu . edu 



J.Y. Yu* 1 

LPSC, Universite Joseph Fourier Grenoble 1, CNRS/IN2P3, INPG, Grenoble, F-38026, France 
& Southern Methodist University, Dallas, TX75275, USA 
E-mail: |yu@lpsc .in2p3.fr 



C. Keppel 

Jefferson Lab, Newport News, VA 23602, USA& Hampton University, VA 23668, USA 



E-mail: keppel@jlab.org 



J.G. Morfin 

Fermilab, Batavia, IL 60510, USA 
E-mail: jnorf in@fnal . gov 



J.F. Owens 

Florida State University, Tallahasse, FL 32306-4350, USA 



E-mail: owens@hep . f su . edu 



T. Stavreva 

LPSC, Universite Joseph Fourier Grenoble 1, CNRS/IN2P3, INPG, Grenoble, F-38026, France 
E-mail: 



stavreva@lpsc . in2p3 . f r 



The neutrino deep inelastic scattering (DIS) data is very interesting for global analyses of proton 
and nuclear parton distribution functions (PDFs) since they provide crucial information on the 
strange quark distribution in the proton and allow for a better flavor decompositon of the PDFs. 
In order to use neutrino DIS data in a global analysis of proton PDFs nuclear effects need to be 
understood. We study these effects with the help of nuclear PDFs extracted from global analyses 
of charged-lepton DIS, Drell-Yan and neutrino DIS data at next-to-leading order in QCD. 
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1. Introduction 

Partem distribution functions are of great importance in contemporary high energy physics. Be- 
cause they encode fundamental information on the structure of hadrons, PDFs are needed for the 
computation of any high energy reaction (at HERA, RHIC, Tevatron, LHC, . . . ) involving hadrons 
in the initial state. 

Many groups perform and regularly update global analyses of PDFs for protons [Q, ^, ^ |]] and 
nuclei [|[ ||, ^]. Although not often emphasized, nuclear effects are present also in the analyses of 
proton PDFs since a number of experimental data is taken on nuclear targets. Most of the nuclear 
targets used in the proton analysis are made of light nuclei in which nuclear effects are expected 
to be small. However, the neutrino DIS data is taken on heavy nuclei such as iron and lead and is 
sensitive to the strange quark content of the proton. It should be noted that a good knowledge of the 
strange quark PDF has a significant influence on the W- and Z boson benchmark processes at LHC. 
Moreover, the neutrino experiments have been used to make precision tests of the Standard Model 
(SM) in the neutrino sector. A prominent example is the extraction of the weak mixing angle sin dw 
in a Paschos-Wolfenstein type analysis. 

Historically, nuclear corrections based on information from the charged-lepton-nucleus (£A) 
DIS data have been applied to the neutrino-nucleus (vA) DIS data. Furthermore, the same correc- 
tion factors have been applied to several different observables in vA-DIS (F2, F3, cross section and 
dimuon production) and at different scales Q 2 . Conversely, it is much more flexible to compute 
nuclear corrections in the Parton Model (PM) using nuclear PDFs taking different observables and 
scales into account. For this reason, we perform global analyses of nuclear PDFs at next-to-leading 
order (NLO) in QCD in a framework closely related to the one used by the CTEQ collaboration in 
Ref. [^]. We determine several sets of nuclear PDFs from a) IA DIS + Drell-Yan (DY) data and 
b) vA DIS data and c) the combined £A DIS + DY + vA DIS data. We analyze and compare the 
resulting nuclear correction factors. 

2. CTEQ Nuclear Parton distribution functions 

The global nuclear PDF framework we use to analyze £A DIS and DY and vA DIS data was intro- 
duced in The parameterizations of the parton distributions in bound protons at the input scale 
of Qo = 1.3 GeV 

*/*(*, Go) =c x Cl (l -xf 2 e c > x (l+e c *x) c s , (2.1) 

where k = u v ,d v ,g,u + d,s,s and 

d(x,Q )/u(x,Qo) = c x c > (1 -x) C2 + (1 +c 3 x)(l -x) C4 , (2.2) 

are a generalization of the parton parameterizations in free protons used in the CTEQ proton anal- 
ysis [|jp. To account for a variety of nuclear targets, we introduce A-dependent fit parameters 

c k ^c k (A) = c kfi + c k ^(l-A- c ^), * = {1,...,5}. (2.3) 

The free proton PDFs are recovered in this framework in the limit A — y 1. From the input distribu- 
tions, we can construct the PDFs for a general (A, Z) -nucleus 

f! A ' Z) (x, Q) = f ff /A M)+ fl' A (x,Q), (2.4) 
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where the distributions of a bound neutron, f^ A {x,Q), are related to those of a proton by isospin 
symmetry. 

We performed a global analysis of IA DIS + DY data within this framework, determining the 
A-dependence of the parameters Ck(A). In this analysis, we applied the same standard kinematic 
cuts Q > 2 GeV and W > 3.5 GeV as in [^] leaving 708 data points and we obtained a fit with 
^ 2 /dof = 0.946 with 32 free parameters (for further details see 

To extract the nuclear effects it is useful to define a nuclear correction factor R which is the 
ratio of the observable (O) using nuclear PDF over free PDF 

OlnucPDF] 

1 J 0[freePDF] v ' 

We compared the nuclear correction factors R for the structure function [F 2 M ] in £A DIS with the 



structure function [F 2 Y ] in vA DIS. As first observed in [11], the nuclear correction factor 7?[F 2 



does not describe the NuTeV vFe data well. This raises the question whether the nuclear correc- 
tions in (A DIS and vA DIS are different. 

For definite conclusions, we set up a global analysis of the combined IA DIS + DY + vA DIS 
data where we included exclusively the neutrino DIS cross-section data coming from the NuTeV 
(iron) and Chorus (lead) experiments, respectively. Here we applied the same kinematic cuts as in 
the first analysis of the charged-lepton data leaving us with 3134 neutrino DIS cross-section data. 
In a fit to only the neutrino data we obtained a # 2 /dof of 1.33 with 34 free parameters (for further 
details see @). 

As expected, the nuclear correction factors using NPDFs extracted from the vA DIS data 
exhibit clear differences with the corresponding nuclear collection factors using NPDFs extracted 
from IA DIS + DY data which are especially marked at low and intermediate Bjorken x. 

Analyzing both data sets in a combined global analysis runs into the problem of an imbalance 
of the number of data points between the two data sets, since the vA DIS data points would be 
dominant compared to the £A DIS and DY data. Therefore, we introduced a weight parameter w to 
combine the data sets 

X 2 = £ xf + £ n 2 . (2.6) 

/ ± A data vA data 

In the case w = 0, only the HA DIS and DY data are included and w = °° uses only the vA DIS 
data. Varying the weight w, we tried to find a compromise fit which would describe both, IA DIS + 
DY data and vA DIS data, reasonably well. We refer to Table II in [[K]] for the resulting % 2 of the 
compromise fits with weights w = 0,1/7,1/2, 1,°° and to Figure 1 in [ |10| ] for the corresponding 
nuclear correction factors. As one can see, with increasing weight w the description of the IA DIS 
data gets worse and the one of the vA DIS data improves. 

In order to judge quantitatively on how well the compromise fits describe the data we use the 
X 2 goodness-of-fit criterion introduced in [ |12"[ ^J. We consider a fit to be a good compromise if its 
X 2 for the £A DIS + DY + vA DIS data is within the 90% confidence level of the fits to a) only £A 
DIS + DY data and b) only vA DIS data. We define the 90% percentile £ 90 used to define the 90% 
confidence level, by 

P(x\N)d X 2 = 0.90, (2.7) 

o 
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where N is the number of degrees of freedom and P{% ,N) is the probability distribution 

( y 2)N/2-l - X 2 /2 

p tf^= %mm ■ (2 ' 8) 

We can assign a 90% confidence level error band to the % 2 of the fits to the £A DIS + DY and to the 
vA DIS data 

x f ±A = 638 + 45.6, XvA = 4192 + 138. (2.9) 



Comparing the results of the fits with different weights, listed in Table II in [10], we conclude that 
there is no good compromise which is compatible with both 90% confidence level given in Eq.(2.9) 
(see details in [|l0|]). 

3. Conclusion 

After performing a thorough global NPDF analysis of the combined charged-lepton DIS, DY and 
neutrino DIS data, we find that nuclear correction factors in vA DIS and IA DIS are different and 
there is no good compromise fit to the combined £A DIS, DY and vA DIS data. This has important 
consequences for global analyses of proton and NPDF, for models explaining the nuclear effects, 
and precision observables in the neutrino sector. The nCTEQ nuclear PDFs described here as well 



as the ones obtained in [ 13] are available at [ 14] 
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